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Blind Adaptive MIMO Receivers for Space-Time Block-Coded DS-CDMA Systems in 
Multipath Channels Using the Constant Modulus Criterion 
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Abstract — We propose blind adaptive multi-input multi-output 
(MIMO) linear receivers for DS-CDMA systems using multiple transmit 
antennas and space-time block codes (STBC) in multipath channels. A 
space-time code-constrained constant modulus (CCM) design criterion 
based on constrained optimization techniques is considered and recursive 
least squares (RLS) adaptive algorithms are developed for estimating 
the parameters of the linear receivers. A blind space-time channel 
estimation method for MIMO DS-CDMA systems with STBC based 
on a subspace approach is also proposed along with an efficient RLS 
algorithm. Simulations for a downlink scenario assess the proposed 
algorithms in several situations against existing methods. 

Index Terms — DS-CDMA systems, MIMO systems, space-time block 
codes, blind adaptive algorithms, interference suppression. 



I. Introduction 

THE ever-increasing demand for performance and capacity in 
wireless networks has motivated the development of numerous 
signal processing and communications techniques for utilizing these 
resources efficiently. Recent results on information theory have shown 
that higher spectral efficiency (T), J2] and diversity (3], can be 
achieved with multiple antennas at both transmitter and receiver. 
Space-time coding (STC) techniques can exploit spatial and temporal 
transmit diversity J3J, |4), 0. The problem of receiver design for DS- 
CDMA systems using multiple transmit antennas and STBC has been 
considered in recent works ll6l- H151 . however, there are still some 
open problems. One key issue is the amount of training required 
by MIMO channels which motivates the use of blind techniques. In 
addition, the existing blind MIMO schemes d, (9), [12), Q3), (B) 
are susceptible to the problem of mismatch resulting from imperfect 
channel knowledge and this calls for a robust approach. Prior work 
on blind techniques for MIMO CDMA systems is limited to sub- 
space |9| and constrained minimum variance (CMV) approaches |8|, 
[12 1, 1 14 1, [15], which are susceptible to the problem of signature 
mismatch. The code-constrained constant modulus (CCM) approach 
has demonstrated increased robustness and better performance than 
CMV techniques 1 161 - 1231 for single-antenna systems although it has 
not been considered for MIMO systems. 

The goal of this work is to propose blind adaptive MIMO receivers 
for DS-CDMA systems using multiple transmit antennas and STBC 
based on the CCM design in multipath channels. In the proposed 
scheme, we exploit the unique structure of the spreading codes and 
STBC to derive efficient blind receivers based on the CCM design 
and develop computationally efficient RLS algorithms. The proposed 
design approach for MIMO receivers requires the knowledge of the 
space-time channel. In order to blindly estimate the channel, we 
present a subspace approach that exploits the STBC structure present 
in the received signal and derive an adaptive RLS type channel 
estimator. We also establish the necessary and sufficient conditions 
for the channel identifiability of the method. The only requirement 
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for the receivers is the knowledge of the signature sequences for the 
desired user. 



II. Space-time DS-CDMA system model 

Consider the downlink of a symbol synchronous QPSK DS-CDMA 
system shown in Fig. 1 with K users, N chips per symbol, Nt 
antennas at the transmitter, N r antennas at the receiver and L p 
propagation paths. For simplicity, we assume that the transmitter 
(Tx) employs only N t — 2 antennas and adopts Alamouti's STBC 
scheme |3|, although other STBC can be used. In this scheme, for 
i = 1, . . . ,P two symbols 6fc (2i — 1) and bk(2i) are transmitted from 
Txl and Tx2, respectively, during the (2i — l)th symbol interval 
and, during the next symbol interval, —b* k (2i) and b k (2i — 1) are 
transmitted from Txl and Tx2, respectively. Each user is assigned 
a unique spreading code for each Tx, which may be constructed 
in different ways (7), JH]. We assume that the receiver (Rx) is 
synchronized with the main path, the delays of the channel paths 
are multiples of the chip rate, the channel is constant during two 
symbol intervals and the spreading codes are repeated from symbol 
to symbol. The received signal at antenna m after chip-pulse matched 
filtering and sampling at chip rate over two consecutive symbols 
yields the M -dimensional received vectors 

K 

r(2» - 1) = Yj Akb k {2i - l)C£h^ + A k b k (2i)C 2 k h 2 m 



+ *Jfc, m (2*-l)+n m (2i-l) 

K 

r(2i) = X>* 6 *( 2i " l)C*hm - Akbl^Clhl 

k=l 

+ V k , m (2i) +n m (2i) 
i = 1 P, m= 1 N r 



(1) 



where M = N+L p — 1, n„ 



[nr(i) ... riM{i)] T is the complex 



Gaussian noise vector with mean zero and £[n m (i)n^(i)] = er 2 I, 
where (.) T and (.) H denote transpose and Hermitian transpose, 
respectively. The quantity E[.\ stands for expected value and the 
amplitude of user k is A k . The channel vector for the users' signals 
transmitted from each transmit antenna n t {nt = 1,2) and received 



at the m-th receive antenna are h„ (i) — [h^ (i) . 



and r/ rn (i) is the intersymbol interference at the mth receive antenna. 
The M x L p convolution matrix CZ* contains one-chip shifted ver- 
sions of the signature sequence for user k and each transmit antenna 
given by s™' = [a"'(l) . . .a^(Af)] T (the reader is referred to Jl], 
1 16|, 1 20 1 for details on the structure of C£'). The received data in (1) 
organized into a single 2M x 1 vector y m (i) = [r T (2i — 1) r T (2i)] T 
within the ith symbol interval at the mth receive antenna is 



ym(») = A ^k{2i - l)C k g m {i) + A k b k (2i)C k g* m {i) + r, k (i) + n(<) 

fe=i 

K 

= Y* k w + * k w + n ^ w + n w 



(2) 
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Fig. 1. Proposed space-time system: schematic of the fcth user (a) Transmitter and (b) Receiver . 



where 



Ck 



a m {i) = 



Cl andC?. eTL Mx 



o , o C; 

o cl ' k -cl 

(3) 
(4) 

The 2 M x 1 received vectors y m (i) are linearly combined with the 
2M x 2 parameter matrix Wfc im (i) of user of the mth antenna at 
the Rx to provide the soft estimates 



where v is a constant to ensure the convexity of ([8} and (|9}, which 
r is detailed in Appendix I along with the convergence properties. The 
P 1 ??%)osed approach is to consider the design problems in (|8] and 
l[9} via the optimization of the two filters Wfe jm (i) and \Vk,m(i) 
l)iril a simultaneous fashion. The optimization of each filter aims to 
) suppress the interference and estimate the symbols transmitted by 
each transmit antenna. The expressions for the filters of the space- 
time CCM linear receiver are derived using the method of Lagrange 
multipliers |26| and are given by 



= W^ m (i)y m (i) = [Zk,m(i) Zk,m(i)Y 



(5) 



By collecting the soft estimates z m (i) at the Rx, the designer can 
also exploit the spatial diversity at the receiver as 



«(o = E 



(6) 



where ot m (i) — diag(« m .i(i), a m ,2{i)) are the gains of the 
combiner at the receiver, which can be equal leading to Equal Gain 
Combining (ECG) or proportional to the channel gains as with 
Maximal Ratio Combining (MRC) (29t 

III. Space-time linearly constrained receivers based on 

THE CCM DESIGN CRITERION 

Consider the 2M -dimensional received vector at the mth receiver 
y m (i), the 2M x2L p constraint matrices Ck and C k that were defined 
in (3) and the 2L P x 1 space-time channel vector g m (i) with the 
multipath components of the unknown channels from Txl and Tx2 
to the mth antenna at the receiver. The space-time linearly constrained 
receiver design according to the CCM criterion corresponds to deter- 
mining an 2M x 2 FIR filter matrix W*. )tn (i) = [vtk,m{i), Wfc jtn (i)l 
composed of two FIR filters Wfc, m (-i) and Wk, m (i) with dimensions 
2M x 1. The filters Wfe, m (i) and Wfc, m (i) provide estimates of the 
desired symbols at the mth antenna of the receiver as given by 

bfc(i)=Bgn(»[w£ m (i)ym(»)]) +j sgn (9f[w£ m (i)y m (i)]) 

(7) 

where sgn(-) is the signum function, selects the real component, 
Ss(.) selects the imaginary component and Wfc, m (i) is designed 
according to the minimization of the following constant modulus 
(CM) cost functions 

JoM(w fcjm (i)) = £[(|w£ m (i)y m (i)| 2 - l) 2 ] (8) 

J C A/(w ft , m (i)) = £[(|w£ m (i)y m (i)| 2 - l) 2 ] (9) 

subject to the set of constraints described by 

Ck w fe , m (i) = v g m (i), Cfw fc , m (i) = v g* m (i) (10) 



w*, m (i + 1) = Rfc,L« [d»,m(») - C k (C%R k X(i)C k ) 1 

■ (Ck Rfc,L(*) d fc,m(«) - v ffrnW)] 

w fc , m (i + 1) = R^ m (i) [d*, m (<) - C k (C k R^(i)C fc ) _1 



i (i)d fe , m (i) - f flfj 



.«)] 



(11) 



(12) 



where R fc , m (i) = E[\z k , m {i)\ ym(i)y m (i)] and R fe , m (i) = 
E{\zk t m(i)\ 2 ym(i)ym(i)] are correlation matrices, where dfc, m (i) = 
#[4,m(*)ym(«)] and d fc , m (i) = £7[^ >m (i)y m (i)] are cross- 
correlation vectors, which are originated from the proposed optimiza- 
tion problem. The expressions (II It and (I12> require matrix inversions 
which lead to a complexity 0{{2M) i ). It should also be remarked 
that i ll It and l l 1 2t are functions of previous values of the filter and 
therefore must be iterated in order to reach a solution. Since i ll It and 
d 1 2b assume the knowledge of the space-time channel parameters, 
channel estimation is required. 



IV. Space-time channel estimation 

In this section, we present a method that exploits the signature se- 
quences of the desired user and the structure of STBC for blind chan- 
nel estimation. Consider the received vector y m (j) at the mth Rx, 
its associated 2M x 2M covariance matrix R m = E[y m (i)y^(i)], 
the space-time 2M x 2L P constraint matrices Ck and Ck given in 
(6) and the space-time channel vector g m (i). From (5) we have that 
the fcth user space-time coded transmitted signals are given by 

x fc (i) = A k bk(2i- l)C k g m {i), x fe (i) = A k b k {2i)C k g* m {i) (13) 

Let us perform singular value decomposition (SVD) on the space- 
time JM x JM covariance matrix R m . 



R„ 



E E[x h (i)x? (t)] + B[xfc(i)*fc (*)] + E[ri h (i)Vk (<)] + ^ 



[V s V n ] 



6T 2 I 



rr 2 I 



[V s V n ] H 



(14) 
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where V s and V n are the signal (that includes the ISI) and 
noise subspaces, respectively. Since the signal and noise subspaces 
are orthogonal [24], [25], we have the conditions V^Xfe(i) = 
V%C k g m (i) = and V^x fe (i) = Y%C k g* m {i) = and 
hence we have Q, = g m (i) H C k H V n V% C k g m (i) = and fi = 
g^(i)C fc V 'nVn C k g* m {i) = 0. From these conditions and taking 
into account the conjugate symmetric properties induced by STBC 
1 8 ], it suffices to consider only Q, which allows the recovery of g m {i) 
as the eigenvector corresponding to the smallest eigenvalue of the 
matrix Cfe H V„V^Cfc, provided V n is known. To avoid the SVD on 
R m and overcome the need for determining the noise subspace rank 
that is necessary to obtain V n , we resort to the following approach. 
Lemma: Consider the SVD on R m as in (14), then we have: 

lim (R m /a 2 )" p = V n V^ (15) 

p — >oc 

Proof: Using the decomposition in (14) and since I + A s /a 2 is a 
diagonal matrix with elements strictly greater than unity, by induction 
we have as p -> oo that (R m /a 2 ) _p = V n V^. 

To blindly estimate the space-time channel of user k at the mth 
antenna of the receiver we propose the optimization: 

Sm(«) = arg mill g"(i)C k H R^ p {i)C k g m (i) (16) 

subject to | \g m (i) 1 1 = 1, where p is an integer, R m (i) is an estimate 
of the covariance matrix R m (i) and whose solution is the eigenvector 
corresponding to the minimum eigenvalue of the JL P x JL P matrix 
Cfe H R rn p (i)Cfc that can be obtained using SVD. The performance 
of the estimator can be improved by increasing p even though our 
studies reveal that it suffices to use powers up to p = 2 to obtain 
a good estimate of V n V^. For the space- time block coded CCM 
receiver design, we employ the matrix Rfc,m(i) instead of R m to 
avoid the estimation of both R m and H. k , m (i), and an equivalence 
of these matrices is established in Appendix II. An analysis of the 
capacity of the system and the necessary and sufficient conditions for 
the method to work is included in Appendix III. 

V. Blind Adaptive RLS Algorithms for Receiver and 
Channel Parameter Estimation 

In this section we present RLS algorithms for estimating the 
parameters of the space-time receiver and channel as described in 
Sections III and IV, respectively. 

A. RLS Algorithm for CCM Receiver Parameter Estimation 

Considering the expressions obtained for w&, m (i) and vf k . m {i) 
in dl 11 1 and d!2t , replacing E[-] with time averages, we can develop 
an RLS algorithm through the recursive estimation of the matrices 
R^Lj T- 1 m (z) = (CfR-; n ( l )C fc )- 1 and fj^o = 

(C k R ; rJ rl (i)Cfe)~ 1 using the matrix inversion lemma (MIL) and 
Kalman RLS recursions |26|. The space-time CCM linear receiver 
estimates are obtained with 

w fe , m (i + 1) = Rfc.L(i) [d fc , m (i) - C k T^ m {i) 

( ~-i - M (17) 

• \C" R fcj ^(i)d fc , m (i) - v g m (i)JJ 

- - _1 r- - - -i 

Wfc, m (i + 1) = Rfe, ro (i) dfe, m (i) - C fc i\ m (i) 

f-n*-i - * M (18) 

■ (C k R feim (j)d fc , m (i) - v g m (i)JJ 

where 

dk, m (i) = ad fc , m (i - 1) + (1 - a)zk tm (i)y m (i) (19) 
dfc, m (i) = ad k>m (i - 1) + (1 - a)zfc, m (i)ym(i) (20) 



correspond to estimates of dk, m (i) and dfc jTrl (irrespectively. In terms 
of computational complexity, the space-time CCM-RLS algorithm 
requires 0((2M) 2 ) to suppress MAI and ISI against 0((2M) ) 
required by i ll It and ( I12t . 

B. RLS Algorithm for Space-Time Channel Estimation 

We develop an RLS algorithm for the estimation of the space- 
time channel g m (i) at the mth receive antenna. The proposed RLS 
algorithm avoids the SVD and the matrix inversion required in 
(19) via the MIL and a variation of the power method used in 
numerical analysis 11271 . Following this approach, we first compute 
the inverse of the matrices R k ] n and R,T m with the MIL, as part 
of the space-time receiver design. Then, we construct the matrices 

r fe , m (i) = Cf R^(i)C fe and f fc , ro (i) = C k A k ^(i)C k . At this 
point, the SVD on the 2L P x 2L P matrices Tk,m(i) and I\ m (i) 
that requires 0(L p ) is avoided and replaced by a single matrix- vector 
multiplication, resulting in the reduction of the corresponding compu- 
tational complexity on one order of magnitude and no performance 
loss. To estimate the channel and avoid the SVD on I\ vm (i) and 
rfe,m(«), we employ the variant of the power method introduced in 

da 

ff m (0 = (i-7(0rfc,m(*))S ra (i-i) (2i) 

where = l/tr\Th >m (i)] and we make g m (i) <- 

g m (i)/||g m (i)|| to normalize the channel. It is worth pointing out 
that due to certain conjugate symmetric properties induced by STBC, 
it is possible to exploit the data record size for estimation purposes 
by using both Tk,m(i) and r\ m (i) and thus the proposed RLS 
algorithm computes 

9 m (i) = (l - 0(O(T*,m(i) + f k, m (i)))g m (i - 1) (22) 

where = l/tr[T ktTrl (i) + I\ m (i)] and the normalization 

procedure remains the same. The algorithm in (22) is adopted since 
it has a faster convergence than (21) due to the use of more data 
samples and spreading codes. 

C. Computational Complexity 

We illustrate the computational complexity of the proposed algo- 
rithms and compare them with existing RLS algorithms and subspace 
techniques, as shown in Table I. The proposed space-time CCM- 
RLS algorithms have a complexity which is quadratic with NtM, 
i.e. the number of transmit antennas N t and proportional to the 
processing gain plus the channel order (M = N + L p ~ 1). The 
complexities of the trained RLS algorithm and the blind CMV-RLS 
|8| are also quadratic with NtM, whereas the complexity of the 
Subspace algorithm of Reynolds et al. |9| is higher due to the 
subspace computation. The complexity of the RLS channel estimation 
algorithm in (21) is 2(NtL p ) 2 + N t L p as it requires the modified 
power method on the sum of the N t L p x N t L p matrices T k>m (i) and 
T k ,m{i)- In contrast, the subspace algorithm of [9] requires an SVD 
or the use of the power method on the NtM x NtM matrix R m . It 
turns out that since M » L p in practice, the proposed algorithm 
is substantially simpler than the one in [9|. 

VI. Simulations 

In this section we evaluate the bit error rate (BER) performance 
of the proposed blind space-time linear receivers based on the 
CCM design (STBC-CCM). We also assess the proposed space- 
time channel estimation method (STBC-CCM-CE) in terms of mean 
squared error (MSE) performance and their corresponding RLS-type 
adaptive algorithms. We compare the proposed algorithms with some 
previously reported techniques, namely, the constrained minimum 
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TABLE I 

Computational complexity of RLS estimation algorithms 

PER SYMBOL. 



N.32, f T=0.0001 , 



Algorithm 


Multiplications 




STBC-Trained 


&(N t M)' z + 2N t M + 2 




STBC-CCM 


5(JV t M) a + 3{N t L P Y 
+4N t M + AN t Lp + 2 


+ 3N?L P M 


STBC-CMV 


4(N t M)' z + 2N'fL p M 
+3N t M + 4N t Lp + 2 


+ 2(N t L p y 


STBC-Subspace 


(N t M) 2 +N?L P M + 
imNtMK + !2N t M - 


3N t M 

f 30K 2 + UK 



variance (CMV) with a single antenna |16] and with STBC |8| and 
the subspace receiver of Wang and Poor without [28] and with STBC 
|9|. The DS-CDMA system employs randomly generated spreading 
sequences of length N = 32, one or two transmit antennas with the 
Alamouti STBC |3| and one or two receive antennas with MRC. The 
downlink channels assume that L p = 6 (upper bound). We use three- 
path channels with powers p]f n given by 0, —3 and —6 dB, where in 
each run and for each transmit antenna and each receive antenna, the 
second path delay (ra) is given by a discrete uniform random variable 
(d. u. r. v.) between 1 and 4 chips and the third path delay is taken 
from a d. u. r. v. between 1 and 5 — T2 chips. The sequence of channel 
coefficients for each transmit antenna nt = 1,2 and each receive 



antenna m = 1,2 



where a lln 



nt 



m l,m 



(i) (I = 0,1,2, ...), 
(i), is obtained with Clarke's model [29|. The phase 
ambiguity of the blind space-time channel estimation method in 1 25 1 
is eliminated in our simulations using the phase of g m (0) as a 
reference to remove the ambiguity and for fading channels we assume 
ideal phase tracking and express the results in terms of the normalized 
Doppler frequency fdT (cycles/symbol). Alternatively, differential 
modulation can be used to account for the phase rotations as in [9| 
or the semi-blind approach of [8| adopted. 

We evaluate the BER convergence performance of the proposed 
RLS algorithms for both receiver and channel parameter estimation 
in a scenario where the system has initially 10 users, the power 
distribution among the interferers follows a log-normal distribution 
with associated standard deviation of 3 dB. After 1500 symbols, 6 
additional users enter the system and the power distribution among 
interferers is loosen to 6 dB. The results shown in Fig. 2 indicate 
that the proposed STBC-CCM receiver design achieves the best 
performance among the analyzed techniques. 

We assess the channel estimation (CE) RLS algorithms with single 
transmit antennas of (28), J25], with STBC of Reynolds et al [9| and 
the proposed space-time RLS channel estimator with STBC given in 
(21) in terms of MSE between the actual and the estimated channels 
using the same dynamic scenario of the first experiment. The results, 
shown in Fig 3, reveal that the proposed space-time channel estimator 
outperforms the single-antenna channel estimator because it exploits 
the information transmitted by 2 antennas. 

The BER performance versus SNR and number of users is shown 
in Fig. 4. We consider data packets of P = 1500 symbols, 2 transmit 
antennas, 1 and 2 receive antennas. We measured the BER after 200 
independent transmissions. A comparison with previously reported 
blind techniques with configurations of 2 transmit antennas and 1 
and 2 receive antennas ((2Tx,lRx) and (2Tx,2Rx)) is shown in Fig. 
4. The curves illustrate that the schemes with multiple antennas at the 
transmitter outperform those with single-antennas and the capacity 
of the system is also increased. With a (2Tx,2Rx) configuration, the 
diversity is further exploited and the proposed STBC-CCM achieves 
a performance close to the MMSE (also with (2Tx,2Rx)) which 
assumes the knowledge of the channel and the noise variance. 
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Fig. 2. BER performance versus number of received symbols in a scenario 
where users enter the system and receivers operate at SNR = E^/No = 15 
dB for the desired user. 
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Fig. 3. MSE performance of channel estimation versus number of received 
symbols in a dynamic scenario where receivers operate at SNR = E^/Nq = 
15 dB for the desired user. 



VII. Conclusions 

We presented blind adaptive space-time block-coded linear re- 
ceivers for DS-CDMA systems in multipath channels. A CCM design 
criterion based on constrained optimization was considered and RLS 
algorithms for parameter estimation were developed. We also derived 
a blind space-time channel estimation scheme along with an efficient 
RLS algorithm. The necessary and sufficient conditions for the 
channel identifiability of the proposed method were established. Sim- 
ulations for a downlink scenario have shown the proposed techniques 
outperform previously reported schemes. 

Appendix 

We present an analysis of the proposed space-time CCM algorithms 
and examine their convergence properties and conditions. Let us 
consider the cost function expressed as (we will drop the time index 
(i) for simplicity) 



Jcm = E 



_(|w^ m y r , 



-1] = J B[| 2fc , m | 4 -2 J B[|2 fe , m | 2 ] + l (23) 
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(a) N-32, K - 10 users, f d T=0.OO01 (b) N=32, SNR = 1 5 dB, f d T=0.0001 

10 ' ■ 10 * t . I I 




SNR (dB) Number of users (K) 

Fig. 4. BER performance versus (a) SNR with K = 10 users and (b) 
number of users (K) at SNR = E b /N = 15 dB . 



Let us recall (2) and define x = YLk=i Xfc = 2fc=i A k b k (2i — 
l)p*:,m- The received data can be expressed by y m — x + x + if + 
n. Since the symbols b k are independent and identically distributed 
(i.i.d.) random variables with zero mean and unit variance, and b k 
and n are statistically independent, we have R = Q + T + cr 2 I, 
where Q = E[xx H ], and T = E[rjr] H ]. Let us consider user 1 and 
antenna 1 as the desired ones, let wi = w and define uu = Akpjfw 
and u = A T P H w = [ui . . . uk\ T , where A = diag(Ai . . . Ak), 
P = [pi, . . . ,Pk], and b = [bi(2i-l) . . .b K (2i- 1)]. Using the 
constraints Ciw = vg and Ciw = vg* , we have for the desired 
user the conditions ui = Aip^w = A\g H and Cfw = vA\g H g. 
Note that the conditions resemble one another and only differ by a 
conjugate term, i.e. Ui =u*. In the absence of noise and neglecting 
ISI, the (user 1) cost function is 

Jcm(w) = E[(u H hh H u) 2 } - 2£[(ubb H u)] + 1 

K K K 

= 8(F + Y, Wkff - 4F 2 J2 Kl* -4F-4J2 KI' 

^ (24) 

where F — mu* = \u-i\ 2 = v 2 A 2 \g H g\. Since u\ = u* and 
\ui\ 2 = \ui\ 2 . The cost functions expressed as above are equivalent. 
Thus, it suffices to examine the properties of only one of them. Let 
us consider the constraint Ciw = vg and rewrite the cost function 
Jcm(w) as 

K 

JcmM = Jcm{u) = S(F + u' H u') 2 -4F 2 ~4Y K| 2 

k=2 

-4F-4(u H 'u') 2 + l 

(25) 

where u' = [u 2 . . . u K f = gw, Q = A' H P' ff , P' = [p 2 . . . p K ] 
and A = diag(A2 . . . Ak)- To evaluate the convexity of v7cm(u'), 
we compute its Hessian (H) using the differentiation rule H = 

a ^H W) a(Jc a A J, (u,)) which yields: 

H = [16(F- i)I + 16u ,jff u'l+16uu' H - 16diag(|-u 2 | 2 . . . K| 2 ) 
1 (26) 
Specifically, H is positive definite if a H Ha > for all nonzero 
a £ C . The second, third and fourth terms of H are pos- 

itive definite matrices, whereas the first term provides the following 



condition for convexity 

vA\\g H g\ > i (27) 

Since u' = Qw is a linear function of w then J7cm(u') 
being a convex function of u' implies that J7ca/(w) = 
Jcm(0w) and Jca/(w) = Jcm(0w) are convex function of w 
and w, respectively, where Q — A' H P' and P' = [pa . . . Pk] ■ 
As the extrema of the cost functions can be considered for small 
noise variance cr 2 a slight perturbation of the noise free case, the cost 
functions will also be convex for small a 2 when vA 2 \g H g\ > j. If 
we assume ideal channel estimation, i.e. = 1, and u = 1, the 

condition will collapse to \Ai\ 2 > 4, which corroborates previous 
results with the constant modulus algorithms [19]. In the case of 
larger values of cr 2 , the designer should adjust v in order to enforce 
the convexity of the cost functions in (8) and (9). 

We discuss here the the suitability of the matrix R.fc, m , that 
arises from the space-time CCM design method, for use in the 
space-time channel estimator. From the analysis in Appendix I for 
the linear receiver, we have for an ideal and asymptotic case that 
(we consider receive antenna 1 and user 1 for simplicity) ut = 
(Aip^)wi ~ 0, for for k = 2, . . . , K . Then, we have that w^y w 
Ai6i(2i-l)wf pi+wf n, and |wf yj 2 w y4 2 |wfpi| 2 +Ai&i(2i- 
l)(wf pi)n fr wi+yli6 1 (2z-l)(pf Wi)wf n+wf nn ff wi. There- 
fore, we have for the desired user : 

R 1 =S[|wfy| 2 yy H ]=/SR + N (28) 

where R = Q + <r 2 I, Q = S[xx H ] = £f =1 |A fe | 2 p fe pf , 
the scalar factor is /3 = Af(|wf^pi| 2 + cr 2 ) and the noise- 
like term is N = A\o 2 [(wf pi (wipf ) + cr (p 1 f wi)(piwf ) + 
([diag(|wi| 2 . . . |w2jv/| 2 ) + ww^] — w H wl)] . 

The conditions presented show that R^ for a general user k can 
be approximated by R multiplied by a scalar factor ft plus a noise- 
like term N, that for sufficient signal-to-noise ratio (SNR) values 
has an insignificant contribution. The same analysis applies to R = 
£/[|w;Fy| 2 yy ], which is given by 

Ri = /3R + N (29) 

\v4iere $ = A\ (|wfpi | 2 + cr 2 ) and N = A\a 2 [(wf pi (wipf) + 
cr(pf wi)(piw 1 T ) + ([diag(|wi| 2 . . . |w 2 m| 2 )+ww h ] -w H wl)]. 
An interesting interpretation of this behavior is the fact that when the 
symbol estimates Zk = w^y and Zk = w^y are reliable and the 
cost functions in (8) and (9) are sufficiently small, then |zfcj 2 and 
\z~k\ 2 have small variations around unity, yielding the approximation 
E[\z k \ 2 yy H ] = E[yy H ] + E[(\z 2 - l)yy H ] <* E[yy H ] = R and 
E[\z k \ 2 yy H ] = E[yy H ] + E[(\z 2 - l)yy H ] Si E[yy H ] = R. 
Therefore, we can employ Rt in lieu of R, since the properties of 
R studied for the proposed space-time channel estimation method 
hold for Rfe. 

We develop an expression of the capacity of the space-time system 
and discuss necessary and sufficient conditions for the identifiability 
of space- time channels and the consistency of the estimates. Let q 3 
and q n denote the signal and the noise subspace ranks, respectively. 
The matrix V^Ch of dimensions r n x JL P will be used for our 
analysis. If the noise subspace V n is the exact subspace then, due 
to V^Cfc^fc = 0, we can verify that the column rank of V^Cfc 
can at most be JL P — 1. In order to have a unique solution (times a 
phase ambiguity) the column rank of V^Cfc must be exactly equal 
to NtL p — 1. Since a column rank of a square matrix is equal to its 
row rank, a necessary condition for a row rank equal to NtL p — 1 
is to have at least NtL p — 1 rows, i.e. q n > NtL p — 1. Since 
q s + q n = NtM this yields 

q s < N T M - N T L P + 1 (30) 
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Consider now the signal subspace rank q s and assume symbol-by- 
symbol estimation in a synchronous downlink system. The number 
of columns of V s (which is an orthonormal basis for Xfe) is q s 
and is composed by the effective spatial signatures of all K users 
transmitted by the Nt antennas, that corresponds to a matrix with size 
NtM x K, and ISI. The ISI corresponds to a matrix with dimensions 
Nt(L p — 1) x K since for our transmit diversity configuration 
we have Nt independent multipath channels with a maximum 
of L p paths each. Assuming that K < NtN, we have that if 
K+2mm{N T (L p -l),K} < NtM— NtL p + 1 = N T (N-1) + 1 
then the necessary condition on q s is satisfied and an upper bound 
for the maximum load of the system with Nt transmit antennas is 

^M«-^- 2 ..,»{(!-^)-.-}] 

(3D 

The above result indicates an increase in the system capacity as 
compared to the result K < N — 2 min{iV/3, L p — 1} for a single 
antenna system reported in 11251 . In order to increase the capacity, a 
designer should choose Nt > 2 . 

Now, we are interested in establishing the sufficient conditions for 
the identifiability of gk- We will follow the approach reported in 
1 24 1 . Let us first consider the effective space-time signatures CkQh 
and Cfe^fe, and rewrite them as 

C k gk = v fc = Xe k /A k , CkQk = v fc = Xe k /A k , (32) 

where = [0, . . . , 0, 1, 0, . . . , 0] T . From the above we have 

fe-l 2M-k 

dim{range(Cfc)nrange(-Y)} = 1, dim { range (Cfc)D range (X)} = 
1 where dim{-} stands for the dimension of a subspace and X = 
Ef=i £[xk(t)Xfc(*)] and X = Ef=i B[xfc(i)x?(t)]. The follow- 
ing theorem establishes the identifiability result for the proposed 
method. 

Theorem: Let C k , X, Ck and X be full-column rank matrices 
and a be an arbitrary scalar. Then the equations 

V^C k g = 0, V^C k g*=0, (33) 

have a nontrivial solution other than agk if and only if the following 
condition holds: 

Condition: There exists the vectors g 7^ agk, and 
such that C k g = X(\ k , C k g = Xq k which are equivalent 
to dim{range(Cfc) D range(Af)} = 1, dim{range(Cfc) PI 
range (X)} = 1 respectively. 

Proof: The equations in d33t can be combined together and written 
as follows V"C k g + V"C h g* = 0. Let the latter condition be 
satisfied with some g ^ agk- As the matrix [X X] is full- 
column rank, the relation ~V„Ckg + V^Ckfj* = directly yields 
[{Ckg) T (Ckg*) T ] £ range(Af X). Hence, from the above we 
have Ckg = X<\ k , C k g = Xq. k - The relations dim{range(Cfe) (~l 
range(A")} = 1 and dim{range(Cfe)nrange(^)} = 1 follows from 
C k g = Xq k , C k g = X<\ k along with d32b and the fact that the 
vectors g k and g are linearly independent. This proves the condition 
in Ckg = X(\ k , Ck§ = Xq k . Now, to prove the sufficiency 
part, let us assume that the condition C k g = Xq k , C k g = X(\ k 
holds. From the latter, we have C k g = Xq k and C k g = Xc[ k and 
[(C k g) T {Ckg*) T ] e range(A" X). It follows that g ^ ag k is a 
solution to V^Cfcg + W^Ckg* = 0. This completes the proof. 
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